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Short-range order (SRO) in Fe-rich Fe-Cr alloys is investigated by means of atomistic off-lattice Monte Carlo 
simulations in the semi-grand canonical ensemble using classical interatomic potentials. The SRO parameter 
defined by Cowley [Phys. Rev. B 77, 669 (1950)] is used to quantify the degree of ordering. In agreement with 
experiments a strong ordering tendency in the Cr distribution at low Cr concentrations (< 5%) is observed, as 
manifested in negative values of the SRO parameters. For intermediate Cr concentrations (5% < cc r < 15%) 
the SRO parameter for the a-phase goes through a minimum, but at the solubility limit the a-phase still displays 
a rather strong SRO. In thermodynamic equilibrium for concentrations within the two-phase region the SRO 
parameter measured over the entire sample therefore comprises the contributions from both the a and a'-phases. 
If both of these contributions are taken into account, it is possible to quantitatively reproduce the experimental 
results and interpret their physical implications. It is thereby shown that the inversion of the SRO observed 
experimentally is due to the formation of stable (supercritical) a'-precipitates. It is not related to the loss of 
SRO in the a-phase or to the presence of unstable (subcritical) Cr precipitates in the a-phase. 
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I. INTRODUCTION 

Iron-chromium steels are of great technological importance 
due to their superior properties at high temperatures and under 
aggressive chemical conditions. In particular, they are used 
in reactor environments because of their good swelling and 
corrosion resistance^ The precise origin of many of these 
beneficial features is still uncertain. Creep resistance is cus- 
tomarily associated with the presence of small a'-precipitates, 
while embrittlement has been related to the presence of larger 
a'-precipitates. The a and a'-phases are the Fe-rich and 
Cr-rich solid solutions into which the body-centered cubic 
(bcc) phase of the alloy decomposes at temperatures below 
about 1000K. At yet higher temperatures a third phase ap- 
pears around equiatomic composition, known as the er-phase. 

The phase diagram for this alloy has been assessed using 
the CALPHAD methodology ^ In this approach the alloy was 
assumed to behave as a subregular solution at temperatures 
below the existence range of the a -phase. In other words, the 
heat of formation is assumed to be positive over the entire 
composition and it is parameterized by a second degree poly- 
nomial in the Redlich-Kister expansion. There is, however, an 
ensemble of experimental observations at relatively low tem- 
peratures in both annealed and irradiated samples, which sug- 
gest a more complex behavior^— 

A breakthrough in the understanding of the microstruc- 
ture of these alloys was made through neutron diffraction 
measurements^ which showed (1) negative Cowley short- 
range order (SRO) parameters (see Eq. (fTJ below) at small 
Cr concentrations indicating strong short-range ordering of 
Cr atoms, and (2) an inversion of the sign of these parame- 
ters with increasing Cr concentration, which seems to imply 
a change of SRO in the solid solution as it approaches the 
solubility limit. The implications of these discoveries for the 
interpretation of creep resistance and embrittlement are appar- 
ent, since SRO is known to affect the mobility of dislocations 9 



and precipitation plays a role in intergranular crackingjAS 

Nonetheless the understanding of the atomistic details of 
the complex behavior of this alloy is still incomplete. This 
situation has motivated a number of recent first-principles 
studies 11-15 which addressed the energetics of this system. 
These calculations revealed a change of sign in the heat of for- 
mation of the solution from the negative to the positive side as 
the Cr concentration increases above approximately 10%.— ~— 
Detailed analysis of the origin of this anomaly showed that, 
while the heat of solution of a Cr impurity in Fe is large and 
negative, magnetic frustration leads to a strong Cr-Cr repul- 
sion, causing the heat of formation to assume large positive 
values as the Cr content increases. 

These advances in the understanding of the energetics of 
Fe-Cr alloys have provided the basis for the development 
of accurate interatomic potential models, which enable large 
scale simulations of the microstructural evolution. Two ap- 
proaches were developed known as the two band model^ and 
the composition dependent modelr^ which both address the 
complex shape of the heat of formation curve at K as deter- 
mined from the aforementioned first-principles calculations. 
A third approach^ to model the finite temperature proper- 
ties of this alloy is based on a cluster expansion fitted to den- 
sity functional theory data which is then used it in combina- 
tion with lattice Monte Carlo simulations to study precipita- 
tion and ordering. Using the composition dependent model, 
we recently studied the implications of the change in sign of 
the heat of formation on the thermodynamic behavior at finite 
temperatures and proposed a modified phase diagram for Fe- 
Cr in the region of low Cr content and temperatures below the 
range of existence of the er-phaseJ2 In fact compared to the 
CALPHAD assessment, 4 which is based on the assumption of 
a standard segregating binary mixture (positive heat of mix- 
ing at all concentrations, compare above), a better agreement 
with the experimental location of the solvus was obtained. 
From our viewpoint three aspects render the study of SRO 
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in Fe-Cr alloys particularly interesting: (1) Experimental 
measurements of SRO parameters of heterogeneous systems 
(such as mixtures of the a and a'-phases) provide only com- 
pound quantities equivalent to averages over the entire sample. 
(2) In Fe-Cr this situation is further complicated by the order- 
ing tendency of the system at low Cr concentrations. (3) The 
recent advances in the development of Fe-Cr interatomic po- 
tentials and the availability of suitable simulation techniques 
provide the possibility to explore the atomistic details of these 
processes. Thereby, we are able to resolve the contributions 
to the experimentally measured quantities and understand the 
microscopic behavior of the material. We discuss the evolu- 
tion of the SRO parameter in a mixture of two phases and 
the evolution of the SRO parameter upon coarsening, induced 
e.g., by aging or irradiation. It is shown that the experimen- 
tally observed inversion of the SRO parameter is due to the 
formation of stable (supercritical) a'-precipitates. It is not re- 
lated to an overall loss of SRO, a loss of SRO in the a-phase, 
or to the presence of unstable (subcritical) Cr precipitates in 
the a-phase. Furthermore, by using the same order parame- 
ter as in the experimental reports, we reconcile an apparent 
discrepancy between experimental data and the model predic- 
tions reported in Ref . [j/7]- 

The paper is organized as follows. Section[fj]introduces the 
order parameters used in this work, while Sect. Hill describes 
the methodology. Section HVI contains the results of the sim- 
ulations organized into three subsections, which address (;') 
the relation between the chemical potential difference and the 
concentration, (;/) the SRO in the a-phase, and (Hi) the evolu- 
tion of the SRO in the two-phase region. In Sect. [V] a simple 
expression for the SRO in the two-phase region is introduced, 
which is subsequently employed to interpret the experimental 
result and to derive the temperature dependence of the SRO 
parameter. SectionlVllsummarizes the conclusions. 



II. MEASURES OF SHORT-RANGE ORDER 

Short-range order in binary alloys is conveniently measured 
by means of the order parameter introduced by Cowley^ 
which has been widely applied in the past (see e.g., Refs. |8|, 
[2lj and H2). The SRO parameter for the fc-th shell of a Cr 
atom is defined as 
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where denotes the number of Fe atoms in the fc-th shell, 

is the total number of atoms in the same shell and cq- 
is the overall (total) Cr concentration. It is instructive to con- 
sider several special cases: (1) If a Cr atom is surrounded ex- 
clusively by other Cr atoms, Z^ = and = 1. (2) If 
one has a completely random distribution the ratio Zp£ /^tot 
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which provides a lower bound to the SRO parameter, in other 
words it indicates the maximum possible degree of short- 
range order at a given composition. 

Previous experimental studies of SRO in Fe-Cr have used 
a SRO parameter specific for the bcc lattice. 8 This parameter 
is defined as a coordination number-weighted average of the 
SRO parameters defined above for the first and second neigh- 
bor shells 
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Other studies^ assumed the SRO parameters a Cr ; and a^ r to 
be the same which leads to similar results. For consistency 
with literature and to simplify comparison we focus on this 
bcc-specific weighted SRO parameter in the following sec- 
tions. In Sect. IIVBI we will however explicitly discuss the 
SRO parameters for different shells, which yields insights into 
the range of Cr-Cr interactions. 



III. METHODOLOGY 
A. Algorithm 

If one is to study the problem of short-range order and the 
formation of precipitates in a segregating alloy while includ- 
ing the effects of local relaxations, one requires an algorithm, 
which is capable of efficiently taking into account both chemi- 
cal mixing and atomic relaxation. The former can be achieved 
by employing a Monte-Carlo (MC) algorithm which allows 
changes in the chemical identity of the particles as will be dis- 
cussed in the following two paragraphs. In principle atomic 
relaxations can be captured by a displacement MC algorithm 
but for the problem at hand molecular dynamics (MD) sim- 
ulations turn out to be more efficient. The basic approach is 
thus to alternate MC moves and MD steps. 

The usual approach to model the effects of chemical mixing 
in heterogeneous solid solutions, is to employ a MC swap al- 
gorithm, in which the global composition is kept constant and 
the system is advanced by swapping pairs of particles of dif- 
ferent chemical identities. The change in total energy is em- 
ployed to decide whether a move is accepted or rejected. For 
single phase systems one can also use a semi-grand canonical 
ensemble (SGC), in which the composition is determined by 
the difference of the chemical potentials A/i of the species, 
and each MC move is equivalent to exchanging particles with 
the external reservoirs of both chemical species. Unlike the 
swap scheme the SGC-MC algorithm can be formulated to 
operate only on one particle at a time. Therefore, the latter 
can be more easily parallelized. 

Within a two-phase region of the phase diagram the rela- 
tion between the concentration c and the chemical potential 
difference has an infinite slope. In principle, the SGC- 
MC algorithm cannot deal with such situations, and in order 
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to keep the global composition within the miscibility gap, one 
must apply an additional constraint. The latter can be formally 
derived as a constraint on the concentration fluctuations, the 
formal derivation of which will be given elsewhere.™ For the 
present study, we implemented the constraint in a very simple 
way, namely as a limitation on the number of particles in the 
reservoirs: When the (externally specified) maximum global 
concentration of solute atoms is reached, MC moves which 
lead to a further increase of the solute concentration are al- 
ways rejected, regardless of their energy, i. e. the maximum 
number of solute atoms is limited. To clarify this approach, 
consider the case in which constraints are imposed on both the 
maximum and the minimum number of solute atoms. In the 
(hypothetical) limit that the difference between the lower and 
upper limit vanishes, this algorithm would become equivalent 
to the swapping algorithm. We recall here that the primary 
motivation behind the development of this algorithm is that it 
can efficiently be parallelized. 

To summarize, for values of A/i within the a-phase, our 
simulations start with pure Fe samples and by alternating MC 
moves with MD steps the sample loads solute atoms and re- 
laxes the lattice in the solid solution until the equilibrium com- 
position is reached and stationary fluctuations reach canonical 
values. If A/i exceeds the stability range of the a-phase, the 
simulation still starts from a pure Fe sample and first estab- 
lishes a solid solution in the a-phase. Since the formation of 
a supercritical nucleus is linked to a nucleation barrier, pre- 
cipitation does not occur immediately but only after some res- 
idence time. The system therefore evolves initially in a sin- 
gle phase. However, as soon as the spatial fluctuations in the 
composition have lead to the formation of supercritical nu- 
clei, they will grow rapidly and the global Cr concentration 
rises sharply. When the target global composition is reached, 
or equivalently, when the volume fraction of a'-precipitates 
reaches the target value, the concentration constraint keeps the 
system fluctuating around the pre-set global composition, and 
coarsening of the precipitates is observed. The residence time 
prior to the formation of stable a'-precipitates depends on A/i, 
i.e. on the chemical driving force for precipitation. If the dif- 
ference of the chemical potentials exceeds the stability range 
of the a-phase only slightly, residence times are rather long 
(in terms of MC steps) and meaningful averages (with regard 
to SRO) can be obtained (see Fig.|4]below for an example). 

As described above, in the early stages of a simulation, 
some a'-nuclei may exist in the a-phase, but they do not 
yet exceed the critical size. Fluctuations in the composition 
eventually lead to the formation of stable (supercritical) a'- 
precipitates, which subsequently coarsen. While this evolu- 
tion cannot be related to a physical time scale, because atoms 
do not move via diffusion but via swaps, it nevertheless rep- 
resents the sequence of steps followed during the real evolu- 
tion of a saturated a-phase. Similarly, SRO develops in the 
a-phase in the course of the simulation, and reaches equilib- 
rium before the a/a'-microstructure emerges, reflecting the 
fact that precipitation is a first order phase transition with a 
nucleation barrier, while order-disorder is not. 



B. Computational details 

We used samples with 40 x 40 x 40 conventional bcc 
unit cells equivalent to 128,000 atoms. Simulations were 
carried out for temperatures between 100 and 900 K and 
chemical potential differences between — 0.45eV/atom and 
0.12eV/atom. To obtain the data presented below, the sam- 
ples were run for a total of 3 x 10 3 MC steps/atom. We have 
found this number to be sufficient to obtain well converged 
results for the SRO parameter. For chemical potential differ- 
ences which are equivalent to Cr concentrations below the sol- 
ubility limit the system equilibrates quickly. Averages are ob- 
tained over all configurations in the Markov chain after the ini- 
tial equilibration period. For values of the chemical potential 
difference which correspond to concentrations slightly above 
the solubility limit the system does not immediately separate 
into two phases due to the nucleation barrier associated with 
this process (compare previous paragraph). Instead one ob- 
serves a transient state during which the system equilibrates 
in the solid solution. Once a supercritical nucleus of the sec- 
ondary phase has formed, it grows very quickly and the sys- 
tem decomposes into two phases. (An example of this behav- 
ior is shown in Fig. [4] below). The existence of this transient 
regime allows us to obtain meaningful averages for concen- 
trations slightly inside the miscibility gap corresponding to 
oversaturated solid solutions. (In Fig.|4]this is the case for the 
first 1000 MC steps). All the data presented in the following 
sections (specifically the data shown in Fig.[2]and Fig. 01 was 
obtained in this manner, i.e. as averages over configurations 
which correspond to the a-phase only. A detailed description 
of a simulation which develops from a single to a two-phase 
system and which illustrates the transient behavior referred to 
above is given in Sect. IIV CI 

The interatomic potential used in this work is taken from 
Ref.Qjjl It is an embedded atom method type potential based 
on the Fe potential by Mendelev et alQ and the Cr poten- 
tial by Wallenius et al.M To model the Fe-Cr alloy, the po- 
tential for the mixed interaction has been constructed to ex- 
actly reproduce the mixing enthalpy curve obtained from first- 
principles calculations^ 

IV. RESULTS 

A. Chemical potential difference vs. concentration 

FigureQ]shows the relation between the equilibrium Cr con- 
centration cci and the difference of the chemical potentials 
A/i = /id- — /LtFe for different temperatures. It is found that 
the dependence cannot be reproduced using the expression for 
an ideal solution, A/i = A/i° + /^Tlnco/ '(1 — cc r ). The 
curves can, however, be reasonably well fit assuming a regu- 
lar solution^ 

A/i = A/i° + 0(1 — 2cq-) + k B T\n CCr (4) 

1 - CCr 

as shown by the lines in Fig. Q] where A/i° and were treated 
as fitting parameters. 
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B. Short-range ordering in the a-phase 

The results for the bcc SRO parameter /3 defined in Eq. 
are shown in Fig. [2] Several observations can be made: At 
small concentrations f3 becomes more negative with increas- 
ing Cr concentration. At low temperature the curves approach 
the theoretical lower limit for the SRO parameter (dashed line 
in Fig. [2j given by Eq. (|2}. For larger concentrations the SRO 
parameter increases, showing a minimum between 8 and 12% 
Cr the location of which is only weakly dependent on tem- 
perature. We conclude that in the entire compositional range 
where only the a-phase exists, its SRO is strong. Furthermore, 
we conclude that the absolute value of the SRO versus com- 
position has a maximum at around 10% Cr, quite independent 
of composition. 

The minimum in the curves occurs due to the onset of a 
saturation of the Cr distribution: Since the interaction energy 
between two Cr atoms in Fe is large and positive* 2 ^ the lowest 
energy configurations maximize the separation between the 
Cr atoms. This is the source of the ordering tendency in this 
system which attempts to maximize the number of non-alike 
neighbors, i.e. zj£' /zf$ in Eq. (fl3 approaches one. 

As the Cr concentration increases, the fraction of configura- 
tions which maximize the Cr-Cr separation decreases. At the 
same time the entropy gained by randomizing the system in- 
creases significantly since the total number of possible config- 
urations increases rapidly. The SRO is therefore both strongly 
temperature and concentration dependent. In the limit of zero 
temperature, i.e. in the absence of entropic effects, the SRO 
parameter approaches its lower bound given by Eq. (|2}. At 
high temperatures it becomes less negative since the entropy 
favors randomization of the system. As the Cr concentration 
increases the entropy gained by randomization eventually ex- 
ceeds the energy gained by ordering causing the minimum in 
the SRO parameter curves. 

This effect is further illustrated in Fig. [3] which shows the 
SRO parameters al^ for the first three neighbor shells. For 
the first shell the calculated SRO parameter at low tempera- 
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FIG. 1. (Color online) Relation between chemical potential differ- 
ence and Cr concentration in a-Fe. The data points were obtained by 
simulation. The lines are fits to Eq. ©■ 
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FIG. 2. (Color online) Dependence of the short-range order parame- 
ter, P, defined by Eq. ifj} on the total Cr concentration. Circles show 
simulation data. Solid lines serve as guide to the eyes. The dashed 
line represents the maximum possible order as given by Eq. (fJJ. The 
solubility limits for each temperature, which is indicated by the open 
circles, were taken from Ref. Qj|. The points which belong to con- 
centrations beyond this limit, i.e. inside the miscibility gap, are indi- 
cated by small filled circles. They have been obtained by averaging 
over configurations of the transient single phase system (oversatu- 
rated solid solution) which is observed prior to nucleation of the ot- 
phase (compare Sect JIIIBl . El: experimental data from Ref.l2ll: E2: 
experimental data from Ref. |22|; E3: experimental data from Ref.@. 



tures follows closely the theoretical limit [Eq. (f2|] correspond- 
ing to a situation in which the first shell around any Cr atom 
is exclusively filled with Fe atoms. With increasing tempera- 
ture the point at which the SRO parameter, , deviates from 
the theoretical limit is shifted towards smaller concentrations. 
For the second and third shells this deviation occurs at yet 
lower temperatures and concentrations. This behavior is cor- 
related with the interaction energy vs. separation curves given 
in Ref. |27j which show a significant drop of the Cr repulsion 
energy with separation. 

As mentioned above the SRO parameter /3 for larger Cr 
concentrations becomes less negative with increasing temper- 
ature. However, for the interaction model used in the present 
work it does not become zero before melting, which would in- 
dicate a completely random solution. It should be pointed out 
that this observation is a property of this model and is not nec- 
essarily related to the behavior of the real alloy, which under- 
goes both magnetic and structural transitions (ferromagnetic- 
paramagnetic, a — 7 — 5) before melting, transitions that are 
not described by the present empirical potential. 

Figure [2] includes several experimental data points. 28 They 
indicate a transition of the SRO parameter through a mini- 
mum and subsequently an increase up to positive values. The 
occurrence of a minimum is seemingly similar to the simula- 
tion results but the two minima have different physical origins: 
The experimentally determined SRO represents an average 
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MC steps/atom 

FIG. 4. (Color online) Evolution of the average concentration and 
SRO parameter during the course of a simulation (T = 500 K, Afx = 
0.06 eV/atom). The points A and B refer to the configurations shown 
in Figures[5]and|6] 
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FIG. 3. (Color online) Short-range order parameters for the first 
three neighbor shells as defined in Eq. Q}. The SRO parameter dis- 
plays a minimum as a function of Cr content and becomes less neg- 
ative with increasing temperature. The latter effect is the more pro- 
nounced the larger the shell radius. The dashed line shows the lower 
bound for the SRO order given by Eq. $2% 



system in equilibrium would be in the a'-phase, whereas the 
constraint keeps the system inside the two-phase region. This 
particular simulation is suitable for illustrating the transition 
from a single to a two-phase system during the course of the 
MC simulation: The initial sample contains Fe atoms only 
and the SRO parameter j3 measured over the entire sample is 
zero. During the very first steps of the simulation the number 
of Cr atoms increases sharply. Already after about 200 MC 
steps/atom the concentration of Cr settles to a value of about 
12% and the SRO parameter j3 adopts a value of about —0.07 
as shown in Fig. [4] In this stage the entire sample is a Fe-rich 
saturated homogeneous solid solution. 

Figure [5] shows histograms of the number distribution of 
SRO parameters obtained by averaging over every Cr atom in 
the sample after 1100 MC steps per atom (point A in Fig. HI 
which is a representative configuration for the Fe-rich ho- 



over the entire sample, eventually including both the a-phase 
and a' -precipitates. In contrast, the simulation data shown 
in Fig. |2] have been explicitly obtained for the a-phase only 
(compare Sect. III1 Bb and the occurrence of a minimum is an 
intrinsic property of this solid solution (as discussed above). 
Thus, in order to explain the relation between the experimen- 
tally measured SRO which involves the entire sample and its 
values in the a and a'-phases, we must first understand the 
contribution of the a' -precipitates to the measured SRO pa- 
rameter. 



C. Bimodal distribution of the SRO parameters in the 
two-phase region 

In order to illustrate the contributions of the a-phase and 
a'-precipitates to the SRO parameter, we consider one par- 
ticular simulation (T = 500 K, A/i = 0.06 eV/atom). At 
this temperature, the chemical potential difference chosen im- 
plies that without constraints on the total concentration, the 
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FIG. 5. (Color online) Number distribution of SRO parameter 
before (A) and after (B) the formation of supercritically sized a'- 
precipitates. The corresponding configurations are shown in Fig. [6] 
The simulation was carried out at 500 K using a chemical potential 
difference of A/i = 0.06eV/atom. 



6 



v t ':.;V; l :---K'; 



4® v. 



E 
re 
re 

O 

cc 

CO 



(a) configuration A: after 1 1 00 MC steps/atom 




20 30 
Concentration of Cr (%) 



■ 





" -; v : ;;> : 



(b) configuration B: after 1500 MC steps/atom 

FIG. 6. (Color online) Atomic configurations before (a) and after 
(b) the formation of supercritical a'-precipitates. Only Cr atoms are 
shown with a color coding based on their SRO parameter: blue/green 
(dark gray) - negative or small positive SRO parameter (Fe-rich solid 
solution, a-phase), yellow/red (light gray) - positive SRO parameter 
(a'-precipitates, a/a'-interface). The simulation was carried out at 
500 K using a chemical potential difference of A/i = 0.06 eV/atom. 



mogeneous solid solution (pure a-phase). There is one pro- 
nounced peak on the negative side equivalent to the SRO of 
the a-phase. Note that the range of SRO parameter also in- 
cludes some positive values. This is related to the saturation 
of the a-phase, which in equilibrium at this temperature has a 
solubility of about 8% Cr— , and to the formation of unstable 
(subcritically sized) a'-nuclei. For further illustration the con- 
figuration is shown in Fig.[6|a) clearly indicating the absence 
of precipitation. 

The system remains in this metastable state up to about 
1300 MC steps per atom. Then it undergoes a rapid transi- 
tion as indicated both by the SRO parameter (average over the 
entire sample) and the concentration. At this point supercrit- 
ically sized a'-precipitates have formed which subsequently 
grow. Without constraints their size would increase until the 
total concentration of Fe reaches the the value corresponding 
to the given chemical potential difference on the Cr-rich side 
of the phase diagram. The final Fe concentration will there- 



FIG. 7. (Color online) Theoretical prediction for the average SRO 
parameter (/?) as given by Eq. lO using Eq. {2]l to express the con- 
centration dependence of the SRO parameter in the a-phase. Ex- 
perimental references as given in the caption of Fig. [2] The set of 
lines between f3 a — 0.8 and f3 a = 1.0 correspond to different val- 
ues of /3 a schematically indicating the effect of coarsening of the 
a'-precipitates. 



fore assume some value lower than the solubility of Fe in Cr. 
In the present simulations, the Cr concentration is however en- 
forced to be about 20% or lower. The distribution of SRO pa- 
rameters shown after 1500 MC steps/atom (point B in Fig.|4| 
in Fig.[5]now clearly displays two distinct peaks which result 
from the a-phase (j3 < 0.2) and the a'-precipitates (j3 « 1), 
respectively (compare Sect.[V|i. 

The resulting SRO parameter averaged over the entire sam- 
ple is about 0.35. Note that the SRO parameter of the Cr atoms 
in the a-phase drops after the occurrence of a'-precipitates. 
This is mostly due to the higher global Cr concentration which 
enters in Eq. <[T|). It is therefore an effect of the concentra- 
tion dependence entering the definition of the SRO parame- 
ter, not an indication for an actual higher degree of order in 
the a-phase. At the same time the concentration of Cr in the 
a-phase returns to the solubility limit. This behavior corre- 
sponds to the curves in Fig. [2] The atomic configuration after 
1500 MC steps/atom is visualized in Fig. Ob) which clearly 
shows the formation of precipitates and the emergence of a 
two-phase mixture. The color scale furthermore illustrates the 
intermediate SRO parameter obtained for atoms at the a/a'- 
interface. 

This example clearly illustrates the transition from a Fe- 
rich saturated solid solution to a two-phase system. It demon- 
strates how in this range of Cr concentrations the local SRO 
parameter can be used to identify atoms as being either part of 
the a or a'-phase (see in particular Fig.|6|b)). Furthermore, 
it elucidates the different contributions which lead to a given 
measurable SRO parameter and how the latter can be utilized 
to detect the formation of a'-precipitates. In the following 
section, these observations will be generalized and combined 
to obtain a general expression for the average SRO parameter 
as a function of Cr concentration. 
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V. DISCUSSION 
A. Interpretation of the SRO inversion 

The results presented above demonstrate that particular care 
must be exercised when interpreting the SRO parameter in the 
case of a two-phase mixture. In the following, we focus on 
the case of an Fe-rich Fe-Cr alloy at low temperature. 

While the SRO parameter for Cr atoms in the a-phase is 
typically small and negative, the SRO parameter for Cr atoms 
in the a'-phase is close to one. Since experimentally the SRO 
is obtained over the entire sample, the measured value is an 
average over all types of local short-range order. This average 
is the weighted sum of the contributions from the a and Ce- 
phases as well as the a/a'-interface. However, for simplicity 
in what follows we subsume the a/a'-interface contribution 
into the contribution from a', i.e. we treat the SRO parame- 
ter of the a'-phase as if it were dependent on the size of the 
precipitates. This assumption greatly simplified the follow- 
ing treatment and the SRO parameter averaged over the entire 
sample can be expressed as follows 

(fi) = Q /3 Q [cg] + (1 - [eg, 1 ] (5) 

where denotes the total concentration of Cr in the system 
and /3 J is the SRO parameter for the Cr atoms belonging to 
the a or the a'-phase, respectively. The fraction of Cr atoms 
in the a-phase, <fi a , can be derived as shown in the appendix. 

According to Eq. © below the solubility limit (j3) equals 
the SRO parameter of the a-phase, the concentration and 
temperature dependence of which have been discussed in 
Sect. UVBl 

We now need to discuss the terms in the equation for ((3) 
which determine the average SRO parameter inside the mis- 
cibility gap: The thermodynamic boundary conditions require 
the Cr concentration in the a-phase to be equal to the solubil- 
ity of Cr in Fe and reciprocally for the a'-phase. Due to these 
constraints the proportions of a and a' depend only on the to- 
tal Cr concentration (compare appendix). Furthermore, since 
the empirical potential^ used in this study underestimates the 
solubility of Fe in Cr, the a'-precipitates contain nearly 100% 
Cr. The solubility also affects the SRO parameter for the a'- 
phase since it correlates with the ratio Zj£ /Ztot^ . However, 
for Fe-rich conditions this can be safely neglected - again - 
because of the small solubility. This conclusion is in line with 
the observations described in Sect. II V Cl 

In Eq. <(5j we have subsumed the contribution of the a/a'- 
interface into (3 a . As shown for example in Figs.|5]and|6|b) 
the SRO parameter for Cr atoms at the interface is signifi- 
cantly smaller than for atoms in the interior. Therefore, both 
/3 a and (j3) depend on the surface to volume ratio of the pre- 
cipitates. The maximum value for (3 a is obtained for the limit 
of a completely coarsened system for which the contribution 
of the a/a'-interface is practically zero (f3 a « 1). At the on- 
set of precipitation the size of the a'-precipitates is, however, 
small and the surface to volume ratio is large. The effective 
f3 a is therefore reduced with respect to the fully decomposed 
system ((3 a ps 0.8). Hence, depending on the kinetics of the 



system f3 a and (/?} can assume a range of values. Both are 
expected to grow gradually as the system coarsens and gets 
closer and closer to equilibrium. This is illustrated in Fig.[7Jas 
described in the following. 

If one uses the theoretical lower limit for (3 a , which as 
shown in Fig. [2] represents the low temperature limit of the 
simulation data for the a-phase, one obtains an entirely ana- 
lytic expression for (0) . The resulting curves for different val- 
ues of /3 Q are shown in Fig. [7J where the solubility of Cr has 
been set to eg 1 = 10% roughly corresponding to the experi- 
mental solubility at this temperature. All curves are negative 
for small Cr concentrations, then pass through zero and even- 
tually approach saturation. The transition from negative to 
positive values is the hallmark of Cr precipitation and rather 
sharp. As discussed above at a given temperature the only 
variable is the SRO parameter for the a'-precipitates which 
depends on their size distribution. The most extremal value 
is f3 a =1 which is achieved if the number of a/a'-interface 
atoms is negligible with respect to the number of atoms in a'- 
precipitates. The resulting curve is shown by the most upper 
line in Fig. [7] For smaller a'-precipitates the SRO parameter 
j3 a is smaller than one. The resulting SRO parameter, (/?}, 
is shown in Fig.[7Jfor several values between f3 a = 0.8 and 
0.95. Upon aging the average SRO parameter at a given con- 
centration and temperature thus will gradually grow from the 
lower to the upper curves in Fig. [7J It is apparent from the 
figure that the effect of a reduced (3 a becomes notable for 
^ > 0.2. 

It should be stressed that the interpretation of the evolution 
of the average SRO parameter (/3) upon aging is solely based 
on a gradual approach to equilibrium, i.e. the coarsening of 
the precipitates. This conclusion is independent of our MC 
simulations. 



B. Comparison with experiments 

On the basis of the foregoing analysis, it is now possible to 
carry out a quantitative comparison with the experimental data 
points included in Figures [2] [7J and [8] Since all experiments 
considered here correspond to temperatures close to 700 K, 
we focus on the SRO curve calculated at this temperature. 

The experimental data point obtained at 5% Cr corresponds 
to a Fe-rich homogeneous solid solution and the simulation 
data for the a-phase (Fig.|2| as well as the predicted average 
SRO parameter (j3) (Fig.|7J are in very good agreement. 

We suggest that at 10% Cr a'-precipitates are formed in 
the experimental samples which drive the SRO parameter up- 
wards leading to an average SRO parameter value close to 
zero. This trend is reproduced by the theoretical curve which 
yields a transition through zero close to the experimental data 
point. The comparison with the calculated curves shows that 
the inversion of the SRO is very steep and is shifted to slightly 
larger concentration with respect to the solubility limit. 

There are several data points at concentrations close to 15% 
which are scattered between (j3) = 0.0 and 0.2: The measure- 
ments in Ref. Hj] (El in Fig. [H) were carried out on a sam- 
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pie aged for 100 hours at 723 K yielding the largest value for 
the SRO parameter amongthe experimental data points. The 
sample measured in Ref. [8] (E3 in Fig. [HJ was subjected to 
isochronal annealing steps at 703 K. In Ref. |22| (E2 in Fig. [8) 
the authors carried out measurements on different samples be- 
fore and after irradiation. All samples were initially quenched 
from 1473 K and subsequently annealed/irradiated at 723 K. 
Irradiation experiments were carried out using both Ar + and 
Fe + ions and exposure times ranging from 3 to 30 min. Within 
the experimental error bars the SRO parameter was zero for 
the initial sample. It was observed to increase with exposure 
time in a series of bombardments with Fe + . The largest SRO 
was obtained after extended irradiation with Ar + ions. 

Collecting the experimental data and using the above analy- 
sis, the following interpretation for the observations at higher 
Cr concentrations emerges: (1) The samples initially obtained 
by quenching from high temperature show a very small SRO 
parameter close to zero. This suggests that during rapid 
quenching from elevated temperatures both the ordering ten- 
dency and more importantly the precipitation and growth of 
a'-precipitates are suppressed (lowest data point from E2). (2) 
On the other hand, the most extensively aged samples show a 
rather large SRO parameter which is close to the values pre- 
dicted by the theoretical analysis presented above (El). This 
suggests that the a'-precipitates in these samples have already 
reached rather large dimensions (diminishing the contribution 
of the a I ol -interface). (3) Samples subjected to energetic ion 
beams showed an increase of the average SRO parameter with 
the intensity of the irradiation, reflecting radiation induced 
precipitation. 



C. Temperature dependence 

The shape of the SRO curves and the transition through zero 
are affected by temperature because of the T-dependence of 
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FIG. 8. (Color online) Theoretical prediction for the average SRO 
parameter based on Eq. (O using j3 a = 0.9 and the simulation data 
from Fig. [2] for the SRO parameter in the a-phase. The temperature 
dependence of the solubility was taken from Ref. [T^ and is indicated 
by the small arrows at the bottom of the figure. 



both the SRO in the a-phase (Fig. |2]i and the solubility of 
Cr in a-FeJ2 If these effects are taken into account one ob- 
tains the curves shown in Fig. [8] The increase of the solu- 
bility with temperature 19 shifts the transition of the SRO pa- 
rameter through zero to larger concentrations, whereas the 
T-dependence of the SRO in the a-phase (Fig. 12) causes a 
slight upward shift of the average SRO. Note that the tran- 
sition through zero is correlated with but not identical to the 
solubility of Cr. 



VI. CONCLUSIONS 

In this paper we have investigated the short-range order 
(SRO) in Fe-rich Fe-Cr alloys using atomistic simulations 
based on an empirical potential description of the alloy. For 
low temperatures and small concentrations the SRO in the a- 
phase is found to be close to the theoretical maximum possi- 
ble SRO. As temperature increases SRO is reduced due to en- 
tropic effects. For somewhat larger Cr concentrations the SRO 
curves go through a minimum, which occurs close to the sol- 
ubility limit. The occurrence of a minimum reflects the com- 
petition between energy induced ordering and entropy-driven 
randomization. It is, however, important to stress that directly 
at the solubility the SRO is still rather strong. 

In the two-phase region the experimentally measured SRO 
parameter is a mixture of both the SRO in the a-phase and the 
a'-precipitates. The MC simulations show that the two con- 
tributions can be clearly separated. When using a weighted 
average of these contributions the average SRO parameter is 
predicted as a function of concentration in good agreement 
with experiment. The MC simulations in conjunction with 
this simple model demonstrate that the inversion of the SRO 
parameter observed experimentally is a result of the forma- 
tion of stable (supercritical) a' -precipitates. It is not related to 
the loss of SRO in the a-phase, since directly at the solubility 
limit there is still an appreciable degree of SRO. The model 
has also been used to investigate the effect of temperature on 
the SRO parameter which shows the temperature dependence 
of the solubility of Cr in Fe to be the dominating factor. 

The results reported here are anticipated to support future 
research in at least two aspects: (1) As indicated in the in- 
troduction the SRO affects the mobility of dislocations. The 
present paper therefore constitutes a basis for a detailed atom- 
istic study of the mobility of dislocations in Fe-Cr alloys fol- 
lowing a strategy similar to Ref. [U (2) The present paper 
clarifies the contributions to the average SRO parameter in a 
two-phase system and introduces a simple predictive model. 
It thereby provides a better understanding of previous exper- 
imental results and provides the basis for a more elaborate 
interpretation of future experiments. 
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Similarly, the number of Fe atoms in the a' -precipitates is 
connected to the equilibrium solubility of Fe in Cr 



sol iua 
C Fe JV Cr 



Combining equations dA.H4A.3l l one obtains 



(A.3) 



Appendix 

The total number of Fe/Cr atoms is simply the sum of the 
number of Fe/Cr atoms in the a-phase and the a'-precipitates 



Using Ng} 



x s F °e N a)- ( A -4) 



'Cr 



and rearranging finally yields 



1 /„.tot sol 

'/Tr x Fe 

l/x s °}-x 



Fe 



(A.5) 



and 



N, 



a ■ 



(A.l) 



In equilibrium the concentration of Cr in the Fe-matrix is 



given the maximum solubility, 



It can be expressed by 



the number of Fe and Cr atoms in the Fe-rich solid solution 
Cq. = N£J (iVp* + Nq t ) which can be rearranged to give 



c Cr JV 



NZ. with 



[/(I 



(A.2) 



This expression gives the fraction of Cr atoms in the a-phase 
as a function of the total Cr concentration within the misci- 
bility gap. For Cr concentrations below Cq. (above cf? ) (j> a is 
one (zero). In summary, (j> a can be written as 



1 

1/x, 



1/3* 



..sol 



sol 



• Cr 



Cr 



Fe 



„tot ^ sol 
c Cr L Fe 



(A.6) 
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